Brane worlds and dark matter 



Shahab Shahidi* and Hamid Reza Sepangr 

Department of Physics, Shahid Beheshti University, G. C, Evin, Tehran (19839) Iran 

September 29, 2010 



Abstract 

Two problems related to dark matter is considered in the context of a brane world model in 
which the confinement of gauge fields on the brane is achieved by invoking a confining potential. 
First, we show that the virial mass discrepancy can be addressed if the conserved geometrical term 
appearing in this model is considered as an energy momentum tensor of an unknown type of matter, 
the so-called X-matter whose equation of state is also obtained. Second, the galaxy rotation curves 
are explained by assuming an anisotropic energy momentum tensor for the X-matter. 

1 Introduction 

Over the past decade, model theories in which our 4D space-time (brane) embedded as a hyper-surface 
in a higher dimensional space-time (bulk) have been proposed and opened a large, relatively unexplored 
area of research. The interest in the study of gravity in spaces with extra dimensions started to grow 
rapidly with the appearance of the proposal put forward by Randall and Sundrum (RS) [1, 2]. This 
was followed by the derivation of the Einstein field equations on the brane through the use of the 
Gauss-Codazzi equations and Israel junction conditions together with the assumption of Z2 symmetry 
[3]. The cosmological consequences of this work, such as the additional quadratic density term in the 
Friedman equations [4, 5] have been widely investigated. 

Although use of the Israel junction conditions together with Z2 symmetry is the widely accepted 
method amongst the workers in the field for confining the gauge fields to the brane and hence obtaining 
the field equation, there has been concerns expressed over their use in that such junction conditions may 
not be unique [6]. Another point to note is that if more than one extra dimension is involved, there is, 
as yet, no way of obtaining viable junction conditions. Eying such concerns, an interesting model which 
avoids the above questions was proposed in [7] . Here, the confinement of the gauge fields to the brane 
is done by postulating the existence of a confining potential V. One then obtains the field equations 
on the brane without the use of any junction conditions or the assumption of Z2 symmetry. The field 
equations obtained in this way contain a term, Q^ v , which is an independently conserved quantity 
and geometrical in nature. This would allow us to attribute it to what is known as the X-matter. 
The dynamics of test particles confined to such a brane- world at the classical and quantum levels was 
studied in [8]. The static, spherically symmetric solutions of this model was used subsequently to 
explain the galaxy rotation curves [9] . Various other cosmological implication of the model were also 
studied and showed to be consistent with observations, such as the late time accelerated expansion of 
the universe which is related to dark energy content of the universe. [10, 11]. 

One of the important observation which necessitates the existence of dark matter is the so-called 
virial mass discrepancy. The total mass of a cluster of galaxies can be estimated in two ways. Knowing 
the motions of its member galaxies, the virial theorem gives one estimate M v . The second estimate 
is obtained by taking the mass of the individual galaxies separately and adding them up to give M. 
Almost without any exception, the virial mass M v is greater than the total baryonic mass M by 20- 
30% [] 2], so that the existence of dark matter is assumed in order to explain this discrepancy. In [13], 
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this effect is expiained within a RS type brane- world together with the fact that the electric part of the 
Weyl tensor can be decomposed irreducibly into dark matter energy density and pressure, previously 
presented in [14]. The virial mass discrepancy is explained in [15] for a brane- world model theory 
discussed above, and also in the context of DGP brane gravity [16]. 

Recent observations show that the tangential velocity of stars moving around the center of a galaxy 
tend to a constant value as we move away from the center of the galaxy. In order to explain such 
rotation curves various suggestions have been made. For example, in [17, 18] Modified Newtonian 
Dynamics is used to explain the stability of circular orbits, while in [19, 20] the authors use conformal 
symmetry in the context of brane- world gravity to explain the constancy of galaxy rotation curves. In 
this paper we present similar conclusions within the context of the brane- world scenario proposed here 
by invoking the notion of conformal Killing symmetry. 

Calculation of the geometrical term obtained in this model needs the specification of the bulk 
geometry. One may then compute the tensor using the Gauss-Codazzi equation. In the case 
that the bulk geometry is unknown, one of the components of the extrinsic curvature tensor is left 
arbitrary, which may then be obtained using the Einstein field equations [21]. In this paper however, 
we use this arbitrariness and consider the geometric tensor as a perfect fluid for which the energy 
density and pressure is obtained via the Einstein field equations. As is shown in the sections below, 
the following two problems related to dark matter can be explained using different geometrical matter. 
The galaxy rotation curves can be accounted for by defining the tensor Q^ v to be anisotropic with 
the equation of state of the form (56). The virial mass discrepancy can be explained by introducing 
the isotropic geometrical matter. This choice however is somewhat expected in that the virial mass 
problem concerns the calculation of the mass of the cluster of galaxies which is a much larger object 
than the galaxy itself. 

The scope of the paper is as follows. In the next section we obtain the virial mass theorem in 
the context of the brane-world model discussed above. In section three the equation of state of the 
X-matter is obtained with the specific choice of the energy momentum tensor of the X-matter. Section 
four concerns the galaxy rotation curves where we show that the constancy of the velocity is obtained 
by the assumption that the brane admits a conformal symmetry. Conclusions are drawn in section 
five. 

2 Virial theorem 

Wc use the Einstein field equations obtained in [15] in the context of brane-world models without the 
Z2 symmetry. Neglecting the global bulk effects, the Einstein equations read 



where t^ u is the energy-momentum tensor of the ordinary matter on the brane and is a conserved 
geometric quantity which can be considered as an energy-momentum of the X-matter. In this section 
we obtain the generalized virial theorem and show that the X-matter mass is linearly related to the 
virial mass. We take the brane spherically symmetric and static metric as 



Gf+v = S-kGt^v — Ag M „ + Q 



(1) 



ds 2 = -e^ r) dt 2 + e v(r) dr 2 + r 2 (d6 2 + sin 2 0)d<p 2 . 



(2) 



For r, 



we take 




(3) 



while for the X-matter we choose a fluid of the form 



Q 



87rGdiag (-p x , p x , p x , p x ) . 



(4) 
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Now, substituting the equations above in equation (1) we obtain the following field equations on the 
branc 

1^1(1 ^ = -87rG(p eff + p x ) - X, (5) 
1 ~'Jr) 



r 
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( 1 + ^ ) - = 8 ^ G (Pe// + Px ) - A, (6) 

(2// - 2!/' - /iVr + 2,/V + / r) = 8ttG(p$ +p x ) — A, (7) 
pi = 2 drPe ff Px r {Pe ff Pe ff > (g) 

Peff +P ( Jff+Px +Px 

where the last equation is obtained from V M (Q MI y + 8hGt^ v ) = and a prime represents derivative 
with respect to r. 

In order to obtain the virial theorem, it is convenient to work in tetrad formalism. Let us then 
define the following frame of orthonormal vectors 

= eU% e« = eU\, ef = r8% ef = rsin^, (9) 

where g^ejf'e^ = rf a ^ b \ and the tetrad indices are indicated by the parenthesis. In tetrad compo- 
nents the 4- velocity u M of a typical galaxy, with the property v^v^ = — 1, is given by 

w («)= w /* e W ] a = 0,1, 2, 3. (10) 

Finally we introduce /(x^ 1 , v<- a >) as the distribution function of galaxies assumed to be identical and 
collisionless point particles. The relativistic Boltzmann equation in tetrad components is now given by 

[22, 23] 

where 7(&)( c ) = e P?°' e ('6) e ( c ) are the Ricci rotation coefficients. For our choice of the metric, the Boltz- 
mann equation becomes 

df .Vgdf * V V df ( V? , vJ+4 \ df V r ( df Of 

Vr ir + e + e — it, \ Ve ir~ + v vir~ 

or r at) r sm dip \ 2 r I ov r r \ ovg av v 

df df 



where 



cot*(^-^)=0, (12) 



W=vt, vW=Vr, v^=v e , v^=v v . (13) 



V 

If we assume that the coordinate dependence of the distribution function is only through r, equation 
(12) reduces to 

Vl dr I 2 r J dv r r \ 8 dvg 9 dv v 

The spherically symmetric nature of the metric requires that the coefficient of cot# be zero. Now, 
multiplying equation (14) by mv r dv, where dv = —dv r dvgdv v is an invariant volume element in the 
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velocity space and m is the mass of the galaxy, and integrating over the velocity space and assuming 
that / vanishes sufficiently rapidly as the velocities tend to ±00, we obtain 

r i \p ^ + \ p ^ + ^ ^ - ' ^ + W> - 2 - °> (15) 

where p is the mass density, and () represents the average value of the quantity it contains. Multiplying 
equation (15) by 4irr 2 and integrating over the cluster of galaxies, we obtain 

'**P [(v* r ) + (vl) + «>] r 2 dr + \ j\^p [<«?) + (vl)} ^-dr = 0. (16) 

The total kinetic energy of the galaxies is given by 

K = j R 2np [(v 2 r ) + (v 2 ) + (vl)] r 2 dr, (17) 



so that equation (16) reduces to 

2K=\j\«r*p[(v*) + (vl)]f r dr. (18) 
In terms of the distribution function we can write the energy momentum tensor t m „ as [22] 

T/iv = / fmv^Vvdv, (19) 



leading to 

Peff=p(vl), P % f =p(v?), V ^) = p(vD= P (vl). (20) 

Now, by adding the minus of equation (5) to equation (6) and adding twice of equation (7) to the 
result we obtain 



pj p!v' p" p! 



= 4:TrGp(v 2 ) +4nG(p x + 3p x ) - X, (21) 



where we have defined (v 2 ) = (v 2 ) + (v 2 ,) + (v 2 ) + (v 2 ). Let us now see what kind of approximations 
we can make in this formalism. First, we assume that p(r) and v{r) are small so that we can neglect the 
quadratic terms in equation (21). Second, we assume that the galaxies have velocities much smaller 
than the speed of light, so (v 2 ) , (v 2 ) , (v 2 ) <C (v 2 ) « 1. These conditions can be applied to test 
particles in stable circular motion around galaxies and to the galactic clusters. Equation (21) is then 
reduced to 

AirGp = 1 1 A ( r V) + A - 4ttG (p x + 3p x ). (22) 

To go any further, we need to specify the equation of the state for the X-matter which we assume to 
have the following form 

Px = a Px > ( 23 ) 
where a is an arbitrary constant. Equation (22) is therefore reduced to 

4^Gp= (rV) +A-47rG(l + 3o)p Jt . (24) 

Multiplying equation (24) by r 2 and integrating from to r we obtain 

\r 2 ^ - GM(r) + ^Ar 3 - GM x (r) = 0, (25) 
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where we have defined 

M(r) = 4tt / pr' 2 dr', (26) 
Jo 

and 

M x (r) =47r(l + 3a) / p x r' 2 dr'. (27) 
Jo 

Finally, multiplying equation (25) by dM ( r > and integrating from to R and using equation (18) we 
obtain 

W + 2K + W x = 0, (28) 

where 

w = j R GJ\£M dM(r)j (2g) 

W x =-[ R ^l d M(r), (30) 
Jo r 

and 

r 2 dM{r), (31) 

is the moment of inertia of the system. Equation (28) is the generalized virial theorem on the brane 
in the presence of the brane cosmological constant and X-mattcr. If W x = 0, equation (28) reduces 
to the virial theorem with a cosmological constant [23]. In what follows we assume that a/-|. In 
order to obtain the ratio of the X-matter mass to the virial mass, we introduce the following radii [23] 

M 2 

Rv = ]» Mil dM {r) ' (32) 
„ f R r 2 dM{r) 

M 2 

R * " j.^' (34) 

where R v is the virial radius and R x is the radius defined by the X-matter. If we define the virial 
mass as 

GM 2 , 

2K =~5 — > ( 35 ) 



Ry 



and use the following relations 



GM 2 GM 
R v 

the generalized virial theorem (28) reduces to 



W = -— — , W x =— I = MR 2 , (36) 

it,, r£„ 



2 

Ml ' ' \ M I \R V I 3G M 



-,V_ 1+ ^Vi (37) 



Since the contribution of a cosmological constant to the mass of the galaxy is several order of magnitude 
smaller than the observed mass, we can neglect it in equation (37). Moreover we can neglect the unity 
in (37) since M v is much greater than M for most galaxies. Therefore, the virial mass in our model is 
given by 



M V ~M X J-^, (38) 

V H x 

showing that the virial mass is proportional to the X-matter mass. 
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2.1 Estimating the virial mass 

In this section we obtain M x as a function of r and consequently see that the solution is linearly 
increasing with r, pointing to a possible explanation of the mass discrepancy in clusters of galaxies. 
We start with Einstein equations (5)-(8) and since most of the baryonic mass in clusters is in the gas 
form, we assume that the effective energy-density and pressure in t^ v is that of a gas and therefore set 

Peff = Pg(r), P ( eff = Piff = Pg( r )- ( 39 ) 

In majority of clusters most of the baryonic mass is in the form of the intra-cluster gas. The 
gas density p g can be fitted with the observational data by using the following radial baryonic mass 
distribution [13, 24] 

2 ~ — 

Pg(r)=Pv + * > ( 4 °) 

where r c is the core radius, and po and (3 are cluster independent constants. A static spherically sym- 
metric system of collisionlcss particles that is in equilibrium, can be described by the Jean's equation 
[12] 

I M + ^ K - < v ) = -P.(r)§, (41) 

where $(r) is the gravitational potential, and a r and crg >v are the mass-weighted velocity dispersions, 
respectively, in the radial and tangential directions. We assume that the gas is isotropically distributed 
inside the cluster, so oy = a$ tip . The gas pressure is related to the velocity dispersion and gas density 
by P g = Pg&r- By assumption that the gravitational field is weak so that it satisfies the usual Poisson 
equation A$ « 47rp tot , where p to t is the energy density including p g and other forms of matter, like 
luminous matter and the X-mattcr, etc., the Jean's equation becomes 

dp g (r) . d* GM t Qt , , .... 

-dr = -f , '^Ar = —^s-P'^ (42) 

where M tot (r) is the total mass inside the radius r. The observed x-ray emission from the hot ionized 
intra-cluster gas is usually interpreted by assuming that the gas is in isothermal equilibrium. So we 
assume that the gas is in equilibrium state having the equation of state [23] 

P 9 (r) = ^Pg(r), (43) 

Lb 1 1 i tj 

where p = 0.61 is the mean atomic weight of the particles in the gas cluster, and m p is the proton's 
mass. Equation (42) then becomes 

Mtodr) = fcsTg r 2 — hip = 3ksT ^ ^ (44) 
prripG dr pm p G r 2 + rj? ' 

where the second equality is obtained using equation (40). Using equations (26) and (27) we can obtain 
another expression for M to t 

= 4?rr 2 pg + 4?r(1 + 3a y px . (45) 

Combining equations (44) and (45) we obtain the following expression for p x 

.M- 

(46) 



l + 3a 



AirprripG (r 2 + r 2 ) 2 \ r 2 



2 \ - — 
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In the regime r ^> r c , equation (46) is reduced to 



1 



l + 3a 



3k B T g f3 
Aitpm v G 



and by equation (27) we obtain 



Sk B T g f3 _ Anpprf ^_ 3p 
pm p G 3(1 — /3) 



(47) 



(48) 



We now consider the limit of large r. For most clusters the parameter (3 has the value /3 > | [24]. 
Because of the smallness of the second term in (48) in this case, we have 

2,k B f3T g 



(irripG 



(49) 



As can be seen, the X-matter mass increases linearly with r and this is similar to the behavior of dark 
matter in the cluster of galaxies. Let us now make some estimates. First we note that k B T g rj 5 KeV for 
most clusters. The virial radius of clusters is usually assumed to be r2oo, indicating the radius for which 
the energy density pi, of the cluster becomes P200 = ^^Punivj where p U niv = 4.6975 x 10~ 30 /i| g/cm 3 
[24]. We can estimate the critical radius of the X-matter by determining the distance in which the 
X-matter density is equal to P200- Using (47) and noting that the last term may be dropped in our 
limit we find 



1.82V/? 



5 KeV 



hJ Mpc. 



We then have 



9.8 x 10 14 (3 



5 KeV 



(50) 



(51) 



For (3 < | the X-matter mass has a maximum at 



Rmax I 



/ 3fc B /3T g 



\ AirGpofirrip 

The maximum value of the X-matter mass is given by 



GAK 



3k B pTg 



,30 



3/3 



Rr, 



(52) 



(53) 



lim p G J \3(l-f3) 

3 X-matter solutions on the brane 

In this section we find the vacuum solutions of the theory and obtain the equation of motion of the 
X-matter. Assuming that t m „ = 0, the effective vacuum equation (1) is reduced to 



Let us assume that Q^ v has the form of an anisotropic perfect fluid 

Q^ = 8irG diag(- Px ,p| ,p±,p 

with an equation of state 



V 



(54) 

(55) 
(56) 
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Taking metric (2), the vacuum field equations on the brane become 



— - — (1 - ri/) - -j = -A - 8itG Px (57) 



^-(l + rfx')-\=-X + S7rGpl, (58) 



e 



-v{r) 



4r 



-(2// - 2i/ - /iVV + 2//V + n' 2 r) = -A + 8ttG^ . (59) 



The equations above contain 4 unknown quantities so that an extra relation would be desirable. Such 
an equation can be provided by assuming that the brane admits a one parameter group of conformal 
motions so that 

= ^Wg^! ( 60 ) 

where <p{r) is the conformal factor and C represents the Lie derivative. Moreover, we assume the 
following general form for the vector field £ 



Use of equation (60) then leads to 



^+2-^ = ^), (62) 



89 dtp 

In order to solve these equations, we see from equation (65) that 



v'e + 1% = >Kr), (63) 

*-k-*4F> (64) 

r^ + ^e=— ' (65) 

V + cot^ + f! = f), (66) 
r dip 2 

sm = — 5T- ( 6? ) 



e 1 = l<f>, (69) 

where we have set the arbitrary separation constant equal to zero (this can be done easily by shifting 
cj) by a constant) and F(ip) is an arbitrary function. From equation (64) we have £° = £°(i) which, 
upon using equation (62), we obtain 

e = f + A, (70) 

where k is an arbitrary separation constant and A is some other constant. Now, using equation (66) 
we obtain 

£ 3 = -cot0F(^)+G(0), (71) 

where G(9) is some arbitrary function. Since ^ is a killing vector field we can set A = without loss 
of generality, so that the form of the conformal killing vector field £ is given by 
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Now, from equations (62) and (63) we can write the metric components in terms of the conformal 
factor 



otir) = C 2 r 2 exp 



-2k 



dr 



(73) 
(74) 



where C and B are arbitrary integration constants. Upon substitution of these equations into the field 
equations (57)-(59) we obtain the following system of equations for 4>{r) , p x and pj- 



B 2 



r 2 B 2 



— A — 8irGp x , 
-A + 8irGp x , 
= -A + 8ttGp^. 



Assuming A = and equating equations (75) and (76) we obtain 



r<j)<j>' + 2<\> 2 

Now, if one expresses r in terms of (j>, one finds 



— B = 0. 



y/2<p 2 - k~4> - B 2 ' 
where ro is an integration constant and f{4>) is defined as 



f(<j>) = cxp 



V8B 2 + k 2 



tanh 1 



-k + 4(f) 



y/8B 2 + k 2 _ 



The energy density and pressure of the X-matter can now be obtained in terms of < 
8irGp x = RirGp" 



1 v/20 2 -k4>-B 2 



B 2 r 2 



B 2 - 3(p 2 + 2k(j>) , 



and 



1 ^2(j) 2 -k<jy- B 2 
B^r 2 M 



(30 2 - 2B 2 - k 2 ) . 



The equation of state of the X-matter can be obtained from equations (81) and (82) 

1 



B 2 r 2 



(k 2 + 2k(/) + B 2 ) 



(75) 
(76) 
(77) 

(78) 

(79) 

(80) 

(81) 

(82) 

(83) 



We see that in the limit r — > oo the X-matter has an equation of state = —p^ 



4 Galaxy rotation curves 

The results of the previous section is in agreement with the observed behavior of the tangential velocity 
of galaxies; we know from observations that the rotational velocities increase almost linearly from the 
center of a galaxy and approaches a constant value of about 200fc™/s as one moves away from the 
center [12]. In this section we consider the tangential velocity of a test particle which moves in a 
circular time-like geodesic orbit. The Lagrangian of the system is given by [25] 
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v tg (km/s) 
200.0 
199.5 
199.0 
198.5 
198.0 
197.5 



1.0 



1.5 



2.0 



r/r 



Fig UI6 ll Variation of vtg as a function of t/tq with B = 1.000000667, and h = 1. 



where r is the afhne parameter along the geodesic. From equation (84) we find that the energy E = e M i, 
the <f component of the angular momentum of the particle l v = r 2 sin 2 6(p, where a dot denotes 
differentiation with respect to r and the total angular momentum of the particle I 2 = I 2 . + sin 6*) 2 
are conserved quantities. The total angular momentum of the particle can be written in terms of the 
solid angle as I 2 = r 4 fl 2 [25]. The equation of the geodesic orbits can then be written as 



where 



V(r) 



r 2 + V(r) = 0, 



; f E 2 e-» - — 



For stable circular orbits, we must have 

r = 0, 



dV_ 

Or 



0, 



d 2 V 
dr 2 



>0, 



(85) 
(86) 

(87) 



so that the potential describes a minimum of the motion. These conditions lead to the following 
expressions for the energy and total angular momentum as [19] 



E 2 = 



2e» 
2 - rjj/ ' 



I 2 = 



r 3 fi' 
2 - r/i' ' 



(88) 



On the other hand the line element (2) can be written in terms of the velocity, measured by an inertial 
observer far from the source as [251 



ds 2 



where 



-dt 2 (1-v 2 



di'l 
~dt 



(89) 



(90) 



In the case of stable circular orbits we have the following expression for the tangential velocity [26] 



"tg 



r 2 e-» 



di'l 
~dt 



Using the expressions for the conserved quantities E and I 2 , we obtain 



' J tg 



s' 1 I 2 
^E 2 ' 



(91) 



(92) 
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Use of equation (88) then leads to 



vl = T ± (93) 



showing that vt g depends only on the time-time component of the metric. Using equation (73), the 
tangential velocity in terms of 4> is given by 



v 



1-* (94) 



Although we cannot express Vt g explicitly in terms of r, we can obtain the value of the tangential 
velocity at infinity from equation (79). We see from this equation that in the case of 4> — ► 4>i,2 where 



k ± Vfc 2 + 8i? 2 

&,2 = ^ , (95) 

we have r — > oo. If we use the upper sign the tangential velocity at infinity is given by 



I 4k 

^ = r- kwmw - (96) 

We only need to fix the ratio of the constants in equation (96) to be consistent with the observational 
data. For example with the choice B/k = 1.000000667 we have v tg0 o = 200km/s. Figure 1 shows a 
plot of Vtg as a function of r. 

Let us obtain the X-matter mass of the galaxy and its dependence on r. After integrating equation 
(57), we obtain 

2(7 

e-» = l-—M x (r), (97) 

where we have defined 

M x (r)=4ir [ Px (r')r' 2 dr' . (98) 



That is equation (27) with a — 0. Now define a boundary R of the galaxy where pb ~ 0. We then have 

e^ = e-" = l-^, (99) 
R 



where 



AI b = 4tt / p(r')r u dr', (100) 

is the total mass of the the galaxy. By equations (74), (94) and (99), we can find the relation between 
constants of the model in terms of the tangential velocity at the boundary 

R 2 i 2GM b 

— = ~^ (101) 
Now, substituting equations (74) and (94) into equation (97) and using equation (101), we obtain 

M x =^r, (102) 

where we have used the fact that the tangential velocity is much smaller than the speed of light. The 
X-matter mass depends linearly on r even after the boundary of the galaxy is reached, and has a value 
of the order of the baryonic mass. This is in agreement with the observational data. 
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5 Conclusions 



In this paper we have considered the problems of virial mass and rotation curves of galaxies in the 
context of a brane-world scenario which uses a confining potential instead of 1i symmetry and the 
Israel junction conditions in order to confine the gauge fields to the brane. We have shown that these 
problems can be adequately addressed in this model by identifying the conserved geometric quantity 
Q^u with a new kind of matter, also known as the X-matter. The virial theorem was then obtained, 
assuming a perfect fluid form for Q^ u and shown to be related linearly to a geometrical mass obtained 
by Qfiv We also assumed an anisotropic form for in order to find the energy-density and the 
pressure of the X-matter and obtained the tangential velocity of a point particle in the galaxy, showing 
that it behaves linearly with respect to r. The X-matter mass of the galaxy was then calculated and 
shown to be of the order of the baryonic mass of the galaxy which extends beyond its boundary linearly 
with r. 
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